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ANALYTICAL DETERMINATION OF COUPLED BENDING-TORSION VIBRATIONS 
OF CANTILEVER BEAMS BY MEANS OF STATION FUNCTIONS^ 

By Alexander Mendelson and Selwyn Gendlbr 


SUMMARY 

yi method based on the concept of Station Functions is pre- 
sented for calculating the modes and the frequencies of non- 
uniform cantilever beams vibrating in torsion, bending, and 
coupled bending-torsion motion. The method combines some 
of the advantages of the Rayleigh-Ritz and Stodola methods, in 
that a continuous loading function for the beam is used, with 
the advantages of the influence-coefficient method, in that the 
continuous loading function is obtained in terms of the dis- 
placements at a finite number of stations along the beam. 

The Station Functions were derived for a number of stations 
ranging, from one to eight. The deflections were obtained in 
terms of the physical properties of the beam and Station Num- 
bers, which are general in nature and which have been tabulated 
for easy reference. Examples were worked out in detail; com- 
parisons were made with exact theoretical results. For a uni- 
form cantilever beam with n stations, the first n modes and 
frequencies were in good agreement with the theoretically exact 
values. The eflect of coupling between bending and torsion 
was shown to reduce the first natural frequency to a value below 
that which it would have if there were no coupling. 

INTRODUCTION 

The failure of turbine and compressor blades due to vibra- 
tions has led to an increased interest in the study of the 
vibrations of these blades and in the determination of the 
natural modes and frequencies. In such theoretical studies, 
it is usually assumed that the compressor or turbine blade 
acts as a cantilever beam. The calculation of the uncoupled 
modes of arbitrarily shaped cantilever beams has been ex- 
tensively investigated (references 1 to 4), but little work has 
as yet been done on calculating the coupled modes of such 
beams. If the geometry of the beam is such that coupling 
exists, the coupled modes are the actual vibrational modes 
that must be calculated. 

Four general methods are currently in use for calculating 
uncoupled modes and frequencies of nonuniform beams. 
These methods are the Kayleigh-Eitz or energy method 
(reference 1), the Stodola method (references 5 and 6), the 
influence-coefficient method (references 4 and 7), and the 
integral-equation method (references 8 and 9) . For each of 
these methods, computational work can usually be carried out 
in several ways. For example, by the use of influence co- 
efficients the modes and frequencies can be determined by 

1 Supersedes NACA TN 218S, “Analytical Determination of Coupled Bending-Torsion V 
Selwyn Qendler, 1960. 


Mykelstad’s iteration procedure (reference 7) or by matrix 
methods (reference 4). 

Any one of these methods can be extended to the calcula- 
tion of coupled bending- torsion modes. The Rayleigh-Ritz 
method usually requires that the uncoupled modes be deter- 
mined before the coupled modes can be computed. In apply- 
ing either the Rayleigh-Ritz or the Stodola method, great 
difficulty is encountered in accurately determining the higher 
modes, because the lower modes must first be “swept out” 
by the use of exact orthogonality conditions (reference 10); 
the process will otherwise always converge back to the 
lowest mode. The same difficulties arc encountered in the 
integral-equation method. 

The influence-coefficient method reduces the problem to 
one having a finite number of degrees of freedom. The beam 
is divided into n intervals and a concentrated loading is as- 
sumed at the center of gravity of each interval. The solution 
of the resultant determinantal equation gives the first n 
modes. The accuracy of the higher modes is, however, very 
poor; only the first third of the modes and the first half of the 
frequencies are obtained within the usual engineering accu- 
racy. Carrying along so many useless modes greatly in- 
creases the labor involved. 

A straightforward accurate method for determining the 
coupled bending- tors ion modes and the frequencies of non- 
uniform cantilever beams, together with applications of this 
method, was developed at the NACA Lewis laboratory dur- 
ing 1949 and is presented herein. This method is based on 
the use of Station Functions as first discussed in refer- 
ence 1 1 . Incorporated in the method are the advantages of 
the continuous-function deflections of the Rayleigh-Ritz 
and Stodola methods together with the advantages of the 
finite number of degrees of freedom of the influence-coefficient 
method. When the method is applied to a uniform beam, 
the first n roots of the resultant determinantal equation are 
amply accurate for engineering purposes. 

The final determinantal equation is solved herein by 
matrix-iteration methods (reference 4). Any other con- 
venient method may, however, be used and no knowledge 
of matrix algebra is needed to carry out the calculations by 
the matrix method. The work can be done by an inexperi- 
enced computer, as the only operations necessary for determ- 
ining each mode are cmnulative multiplication and division. 
In addition, for the case in which the coupling coefficient 
remains constant along the beam, a simple quadratic 
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formula and a series of curves are presented for determining 
the first coupled mode in terms of the uncoupled modes. 
Examples are developed in detail and comparisons with 
exact theoretical results are included. 

THEORY 

In the usual influence-coefficient methods for solving 
dynamical problems, a continuous body having an infinite 
number of degrees of freedom is replaced by a body having a 
finite number of degrees of freedom. Two principal assiunp- 
tions are then made that introduce inaccuracies into the 
solutions, particularly in the higher modes: (1) The resultant 
of the inertia loads of all the infinitesimal masses in a finite 
interval passes through the center of gravity of that interval; 
and (2) a concentrated load that is the resultant of a dis- 
tributed load produces the same deflection as the distributed 
load. An attempt has been made to reduee the error due to 
the second of these assumptions by the use of weighting 
matrices (reference 12). Although the accuracy is thereby 
increased, the effect of the first assumption is still great 
enough to introduce serious errors (reference 11). 

In order to eliminate these assumptions, Rauscher (ref- 
erence 11) introduced the concept of Station Functions. 
Instead of assuming the inertia loads to be concentrated at 
the centers of gravity of the intervals, the inertia loads 
and, consequently, the deflections are assumed to be con- 
tinuous functions along the beam. The values of these 
continuous deflection functions at the reference stations must 
equal the deflections of the reference stations. The loading 
on the beam is therefore a continuous function of the de- 
flections of the reference stations. Inasmuch as the deflec- 
tions of the reference stations can be computed from the 
loading on the beam, which in turn is available from the 
deflections, the deflections are therefore obtained as functions 
of themselves. This procedure gives n homogeneous equa- 
tions in the n deflections of the reference stations. The 
resultant determinantal equation has n roots for the fre- 
quency; it will be shown that for a uniform beam all these 
roots are sufficiently accurate for engineering purposes if 
the deflection functions are properly chosen. (For coupled 
bending-torsion vibrations, 2n homogeneous equations and 
2n roots are obtained for n stations.) 

The deflection fimctions used must satisfy the boundary 
conditions of the problem and also the condition that, at 
any reference station, the value of the function must equal 
the deflection of the reference station. Although it is always 
possible to find directly a single function that will satisfy 
these conditions, it is more convenient to obtain different 
component functions at each station and to add all these 
component functions together to give the complete deflec- 
tion function. Rauscher (reference 11) calls these compon- 
ent deflection functions Station Functions. For example, 
the complete torsional deflection function for the beam will 
have the following form: 

(«s)=S 

1=1 

where 

z dimensionless distance along beam 

e{z) torsional deflection at distance s from root 

6] torsional deflection at station 


Jj{z) Station Function in torsion associated with station 
(All symbols are defined in appendix A.) 

Each Station Function must satisfy the boundary condi- 
tions of the problem and the following additional conditions : 
(1) At the reference station with which it is associated, the 
Station Function equals the deflection of that reference sta- 
tion; and (2) at all other reference stations, the Station 
Function equals zero. The sum of all these Station Func- 
tions will then give the complete deflection function for the 
beam. The Station Functions and corresponding loading 
functions are derived in appendix B for torsional vibrations, 



bending vibrations, and coupled bending-torsion vibrations of 
an arbitrary cantilever beam. 

Torsional vibrations. — It is shown in appendix B that the 
torsional deflections of the reference stations for a beam 
divided into n intervals of length 6, as shown in figure 1, 
are given by the following system of equations : 

= atjdj ( 1 ) 

<^0 j=i 

where 

g.<= 7^ r IirNik—jk — (2) 

k=i Ok L r=k+i J 

i and j= 1, 2, ... n 

CO frequency of vibration 
S length of interval 

lo mass moment of inertia per unit length about elastic 
axis at root section 

ratio of average mass moment of inertia per unit length 
of interval to mass moment of inertia per unit length 
at root section 

Co torsional stiffness of root section 

Ck ratio of average torsional stiffness of fc*** interval to tor- 
sional stiffness at root section 
The Station Numbers Njk and Mjk are functions only of the 
integers k, j, and n and are defined as 



where fj{z) represents the Station Functions derived in 
appendix B and is given by 

fj{z) = aijZ-\-a2iZ^-}- ■ . . +a(K+i)jS (4) 

The coefficients at, are determined in appendix B by 
satisfying the conditions on the Station Functions. The 
integrals in equations (3) are thus seen to be integrals of 


I I I I I nil 


nil I nil I II Hill II I II nil 


II 


III II 
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simple polynomials and the hmits of integration are integers. 
The Station Numbers and Mjic are therefore rational 
numbers, fimctions only of the integers n, k, and j. These 
numbers have been evaluated and are listed in tables I to 
VIII. 

If the physical properties of the beam under consideration 
are known for each of the n intervals, and will be 
known. The Station Numbers Nj* and can be obtained 
from tables I to VIII. From equation (2), at] can then be 
easily calculated. 

Equation (1) actually represents n homogeneous equations 

• 1 o 

m the n unknown deflections 0(. With — these equa- 

-/qw 

tions can be written as follows: 


(5) 


(ail — ^)^i+ai202'hai303-[- 

a21®l+(a22 — X)02+ a2303-l- 

“31^1 + a3202+(a33 X)03-|- 

• • • +ai„0„ = O 

• • • Ol2n0n — 0 

■ ■ ■ +a3„0„ = O 

amSi-f a„2^2 + a„303-t- . . . 

+ (ann — X)S„ = 0 


For a nontrivial solution, the determinant of the coefficients 
must vanish and the characteristic equation becomes 


ail- 

-X ai2 

ai3 

• . <Xln 

a2i 

a22 — 

X a23 

■ ■ Ol2n 

«31 

«32 

a33— X . 

■ ■ ag,, 

a„i 

a„2 

a„3 • ■ 

• n 


= 0 


( 6 ) 


or 

|X7 — [<Xij]\ — 0 (6a) 

where / is the identity matrix, and [«y] is the dynamical 
matrix. 

Equation (6) can be solved for the n values of X by any 
method available. The method used herein was to obtain 
the values of X as the latent roots of the matrix [ay], which 
is actually the dynamical matrix for the problem. The mode 
shapes are obtained at the same time. 

Bending vibrations. — The bending deflections for the beam 
shown in figure 1 are given by the following system of equa- 
tions (appendix B) : 


-^0 i=l 


(7) 


where 

i and i=l, 2, ... n 

mo mass per imit length of beam at root section 
ms ratio of average mass per unit length of 7'*' interval to 
mass per unit length at root section 
Bo bending stiffness at root section 

Bt ratio of average bending stiffness of k'’’^ interval to bend- 
ing stiffness at root section 


(9) 


The Station Numbers M' N'j^, P'j^, and Q'j^ are func- 
tions only of the integers k, j, and n and are defined by 

(* — 1)"J 9j{2)d2 

J*^ {k — l)^z+^ (A: — 1)®J 

r g}{z)dz 
Jit-1 

N'ik= r zgi{z)dz 
J k-l 

The Station Functions gj{z) are derived in appendix B and 
are given by 

gj{z) = b2jz^+b3jz^+bijz^+ . . . +b^n+3)JZ^"■+^'> (10) 

The integrals in equations (9) are thus seen to be integrals 
of simple polynomials. The Station Numbers M'jk, N'j^, 
P'jk, and Q'jt: are rational numbers, functions only of the 
integers j, k, and n. These numbers have been evaluated 
and arc listed in tables I to VIII. 

If the physical properties of the beam are known for each 
of the n intervals, mt, and Bi, will be known. The Station 
Numbers N' P'^^, and are obtained from tables I 
to VIII; {iij can then easily bo calculated by equation (8). 
The determinantal equation is: 


dll X dl2 

dl3 

• Pin 



d21 

d22 

X d23 

■ Pzn 



dsi 

ds2 

d33 X . 


=0 

(11) 

ditl 

dn2 

dn3 

• ^nn ^ 




or 


where 


|X7-[|3y]| = 0 


(11a) 


mo S'* Cl)''* 

The dynamical matrix is [/3y] . 

Coupled bending-torsion vibration. — The torsional and 
bending deflections due to coupled bending-torsion vibra- 
tions of a cantilever beam are given by (appendix B) : 


6i — w 


Bo 




(12) 


where 


fa 


F=— — 

~ Co mo 

ro absolute magnitude of projection of distance from elastic 
axis to center of gravity on perpendicular to bending 
direction for root section 

r*o radius of gyration about elastic axis at root section 


I I I I II 


I Him nil 


II II 1 1 III iiii I 


nil II 
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The quantities atj and are defined by equations (2) 
and ( 8 ). The quantities yu and are given by 


t=l Oft L T=k+1 J 

hj - ± ^ (^S,{iFj,~ Qjft) i: ^ s. I(i - fc +l)Njr + 


y ( 13 ) 


where 

(3, ft ^ fi, [i-l j 

and fSft is the ratio of the average static mass unbalance of 
the interval to the static mass unbalance at the root 
section. 

The Station Numbers Py* and Qji, are listed in tables I to 
VIII with the other Station Numbers. The determinantal 
equation becomes 


ran— 

X Fai2 

. rai,j 

eryii 

«F7i2 . 

■ «r7i„ 


of stations. 


r a2i 

r Ci22 — ^ • 

. ra2)i 

€F721 

eF722 • 

• «F72„ 


Let 71=1 ; 

aii=A7’ii 








From table 1, Nu 

= 5/12, 

Ta„i 

r««2 

• • 1 cKk n 

X er7„i eT7„2 • 

• er77iB 



a,i = 5/12 





=0 


5ii 

5i2 

. . 

Pn- 

X /3i2 

• • Pin 


and 


521 

^22 

• 52)1 

Pn 

Xu 

to 

1 

■ l^2n 


or 


5nl 

5)12 

. 5„,1 

Pnl 

Pn2 

• • Pnn X 



0 

0 

II 

=ok 
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3 


or 


hij]l = o 


(14) 

(14a) 


where [77 jy] is the dynamical matrix and I is the identity 
matrix. 

The first ?i roots of equation (14) will give the first n 
coupled frequencies. 

APPLICATIONS AND RESULTS 

In applying the previously discussed method, it is necessary 
to determine for a given beam the elements ay, /3y, 7 y, and 
Bij of the dynamical matrices. These quantities will depend 
on the physical properties of the beam and on the number of 
stations chosen. If the physical properties of the beam are 
known, the quantities ay, /Sy, ry, and 5y can be directly cal- 
culated from equations ( 2 ), ( 8 ), and (13). The numbers 
Mj,:, Nit, Pm Qm M'm N'm P'm and Q'jt appearing in these 
equations depend on the number of stations n that are used 
and can be read dhectN from tables I to VIII for any given 
number of stations up to eight. Once these quantities have 
been calculated, equations ( 6 ), (11), or (14) can be solved 
for the frequencies by any method desired. The matrix- 


iteration method used herein is simple and rapid and re- 
quires no particular computing skill. As will be indicated, 
however, the accuracy of equations ( 6 ), (11), and (14) is 
such that relatively few stations need be used, in which 
case it may be convenient to expand the determinants and 
to solve the resultant low-order algebraic equation. 

In order to illustrate the accuracy, this method was applied 
to torsional vibrations, bending vibrations, and coupled vi- 
brations of a uniform cantilever beam. The exact theoretical 
values for torsional vibrations and bending vibrations of 
uniform cantilevers are well known. The exact theoretical 
values for the coupled bending-torsion vibration of a uniform 
beam were calculated (appendix D). A comparison was 
then made between the values obtained by the method 
presented and the exact theoretical values. The number of 
stations used was 1 , 2 , and 3 (n=], n=2, and n=3). The 
comparisons are summarized in table IX. 

Torsional vibration. — For the case of a uniform beam, 
Ck=Tk=^ and equation ( 2 ) becomes 


S r Njk- 


k=l L 


■ik- 


■ 1 ) 

r=t+l 


(15) 


The values of Njk and Mys are given in tables I to VIII. 
The table to be used depends on the choice of the number 


j=1.549 


I Co 

\ I ol^ 


The exact theoretical value for the first torsional frequency 
is 

The percentage error is — 1.4 when only one station is used. 

The mode shape obtained by the method of Station Func- 
tions agrees well with the theoretical mode shape, as is shown 
in figure 2 (a) . 

Let 71=2; then by equation (15) and table II, 

AT I Tir 8.5 57 

a„-Nu+M,2-j^+j2-^ 


0 C 12 — Njid- M 22 — ■ 


31 29 57 

'240'^48~120 


AT , AT 8,8 16 


“22 = N2l-f Af2S 


31 239 _13 

'240'*“240 15 
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(a) First mode, 

(b) First mode, 7i = 2. 

(c) Second mode, ?i = 2. 

FinuRE 2.— Comparison of theoretical mode shapes with mode shapes 



(d) First mode, ?i=3. 

(e) Second mode, ;?,=3. 

(f) Third mode, n=‘i. 


by taking n stations along the beam for torsional vibrations. 


The determinantal equation then becomes 


Let n=3; then by equation (15) and table III, 


which gives 


Therefore 


60‘ 

16 

15 


57 

120 


15 ^ 


= 0 


Xi=1.6214 

X2=0.1953 


»,=4.526^A 


The exact theoretical values are 

-,= I.671^g. 


The precentage errors of the first two modes, for only two 
stations, are found to be 0 and —4. 

The mode shapes are shown in figures 2 (b) and 2 (c). 
Agreement of the first mode with the exact theoretical shape 
is excellent; the second mode agrees fairly well. 


QJii — iVn “h Afi 2 “h A/x 3 — 0.945833 
a 1 2 = A ^2 1 “t" A /22 “I” Af23 = 0 . 9 5 8 3 3 3 
^13 ~ iVsi -j“ A/ 32 -l-A/ 33 = 0 . 520834 
« 21 = A^U + A^i 2 + 2 A/,3 = 1 .033333 
1 .883333 

<X23 — ^31~\~N22'\~^M33= 1 .011 1 13 
0!31 = A^ll-l“fVi2“f'fVl3=^ 1 .012500 
®32 — A^22 -h A^23 = 2 . 02 5000 
0 : 33 = fV 3 i A732-1- fV33= 1 .387 501 


The determinantal equation is 


0.945833-X 0.958333 
1.033333 1.883333 — X 

1.012500 2.025000 


0.520834 

1.011113 

1.387501-X 


= 0 


The solutions are 

Xi=3.6474 
X2= 0.4093 
Xs=0.1599 


Therefore 




0 ) 3 = 7.502 



* 
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(a) First mode, 71=1. 

(b) First mode, n~2. 

(c) Second mode, ti=2. 


(d) First mode, ti=3. 

(e) Second mode, n=3. 

(f) Third mode, 7^=3. 


Figure 3.— Comparison of theoretical mode shapes with mode shapes obtained by taking n stations along the beam for bonding vibrations. 


The exact theoretical values are 

The percentage errors of the first three modes, calculated 
by use of three stations, are found to be 0, —0.5, and —4.5, 
respectively. 

The mode shapes are shown in figures 2 (d) to 2 (f). The 
first two modes agree very well with the theoretical shapes; 
agreement of the third mode is fair. 

This procedure can be carried out as shown for any number 
of stations desired. 

Bending vibrations. — For a uniform beam, 5*=m,,= l and 
equation (8) becomes 

( 1=1 (. r=t+l L\ 

I (1 6) 

Let n= 1 ; 

.■.Pn = P'n-Q'n 


and from table I 

71 31 _ 59 

^""630 1008 720 

Therefore, from equation (7), 

The exact theoretical value is 

The precentage error for just one station is found to be 
-0.65. 

The mode shape is shown in figure 3 (a) and is seen to 
agree very well with the theoretically exact shape. 

Let 71=2; then by equation (16) and table II, 

du=-P'n-Q'„+iiV',2-iM' 13=0.422745 

/3i,=P'3,-Q'3i+iiV'33-|M'32=0.295925 

fe= 2F'u + 2 P'i 3- Q'n - Q'i2+|iV'i3- 1 M'i3= 1 . 1 45 1 6 7 

fe=2F'2, + 2F'22-Q'21-e'22 + |iV'32-|M'22=0.905530 
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The characteristic equation is 


0.422745— X 
1.145167 


0.295925 

0.905530-X 


= 0 


The roots are 


Xi=1.2943 


X 2 = 0.0339 




= 3.516 


co2=21.71 

The exact theoretical values are 
ui = 3.516 


/ B, 

\mol^ 




V 


The percentage errors for two stations are therefore found 
to be 0 for the first mode and — 1.5 for the second mode. 
The mode shapes are plotted in figures 3 (b) and 3 (c). The 
first mode agrees excellently with the theoretically exact 
shape; the second mode agrees fairly well. 

Let n = 3; then by equation (16) and table III, 


|3ii — P' u 2 •^^12+ 2 


4 \ M'i3=0. 270604 

6 6 


The characteristic equation is 
0.270604-X 1.009943 

0.648170 3.266250-X 

0.985135 5.822852 


0.487441 
1.689891 
3.204301— X 


= 0 


The roots are 


Therefore 


Xi=6.5521 

X2=0.1667 

X3=0.0223 


toi = 3.5 1 6 

/ Po 


C 02 22.04 

\mol^ 

co3 = 60.20 

rK 

"V 


The exact values are 


^12 — P '21 9 •^^22+ 2 '^'' 23 ' 0 ^'i' 22 — 0 — 1.009943 

= + + M' 32 -^ M'33= 0.487441 

021 = 2 P''ii-l- 2 P ^12 Q' w ^^2 4■ 

|iV ',2 + 2Ar'i3-| M', 2 -| M',3 = 0.6481 70 

d 22 = 2 P '21 + 2 P' 22 — Q' 21 Q' 22~\~ 

|iV' 22 + 2Af'23-| M' 22 -| M'23 = 3.266250 

^23= 2P%l-t- 2P^32 <2^31 Q's2~1~ 

|A^'32 + 2Af'33-|M'32-|M'33=1.689891 

/33I = 3P^H-1- 3P^12“1“ 3P^13 Q\i Q'i2 (2^13 4" 

|fV',2+4iV'.3-^M',2-^M'>3=0.985135 

^32= 3P'21+ 3P^22“1"3P^23 <2^21 — Q' 22 Q^23 + 

I A’"' 22 + 4 Af '23 - 1 M '22 - ^ M '23 = 5 . 8 2 2 8 5 2 

^ 33 = 3P^31 "1“3P^32 4“ 3 P^ 33 — Q' 31 — Q' 32 — Q^ 33 “t“ 

|7V'32+4iV'33-jM'32-^M'33 = 3.204301 


The percentage errors for three stations are found to be 0 , 
0, and —2.4, respectively. The modes are plotted in figures 
3 (d) to 3 (f). The first two modes are seen to agree very 
well with the theoretical mode shape; agreement of the 
third mode is fair. 


Coupled bending-torsion vibrations. — A uniform beam 
with the following constants was chosen: 


7= — 2=38.56 
^6 


«= 0.8 




193.2 


eP = 


71“ 

^241.5 


The values of Uij and are obtained as previously and are 
the same as given before for n=l, ?i=2, and ti=3. Also, 
because /Sa=Pj;= C*= TOt=P= 1 , equations (13) become 


T/j 


-ik- 


M' jr 

r=A+l 


k^l L 

5„= ^ [(i-* +i) A^,,-b 


c 


¥-(Jc-lf 2k-l 


i)Mj| 


928716—51 2 
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Let 71=1; then the determinant is 


Fan — ^ «F7n 
^li fin ^ 


0.002156-X 

0.111111 


0.001196 

0.081944-X 


=0 


The roots are 


Xi = 0.0837 
X2 = 0.0005 

o>, = 3 . 46 ^^“,- 
wo = 44.7 -x/^ 


The procedure for calculating the e.xact theoretical values is 
derived in appendix D. The exact values are 



The percentage error for the first mode, calculated by use of 
one station, is — 0.9. 

Let n = 2; then the determinant is 


r an — X 

F «12 

eFti, 

eF7i2 

r Qfoi 

Va22 — '>^ 

eF72i 

eF722 

Sn 

^12 

/3.1-X 

fil2 


5o2 

fi-2X 

fi-22 — X 


Substituting the known values and solving for X give for the 
first two roots 

X, = 1.3197 
X,= 0.0412 

and the frequencies become 

The percentage errors for two stations are — 0.3 for the first 
mode and —4.4 for the second mode. 


This procedure can be carried out for any number of sta- 
tions desired. For three stations, the frequencies obtained 
are 



The precentage errors are — 0.3 for the first mode, 0 for the 
second mode, and —1.8 for the third mode. 

The results obtained by the method presented are seen 
to agree very well with the exact theoretical values. 

These results are summarized in table IX, where a com- 
parison is made with the results obtained for uncoupled 
bending and torsional vibrations by use of influence coeffi- 
cients with weighted matrices (reference 12). The values 
using weighted matrices were taken from table I of refer- 
ence 12. It can be seen that for a given number of stations, 
the results obtained by the method presented herein are con- 
siderably better than those obtained by using influence co- 



FitH’RE 4.— Varialion of frequency ratio with coupling: cocfficionl e for several values of 
uncoupled IVtHiuency ratio "V. 

efficients with weighted matrices. In general, it is indicated 
that for a uniform cantilever beam using stations along 
the beam, the first 71—1 frequencies and modes arc in ex- 
cellent agreement with exact theoretical values and even 
the mode is given within the accuracy with which the 
physical properties of the material are known. For a tapered 
beam, more stations may be required, depending on the 
amount of taper. The number of stations required to give 
satisfactory accuracy is listed in table X. A comparison is 
made by using weighted influence coefficients; the values 
are taken from table II of reference 12. 

The first vibrational frequenc}^ is given approximately by 
equation (C2) (appendix C) when coupling exists between 
bending and torsion; it is plotted in figure 4. In order to 
check these curves, the exact solution was obtained (appen- 
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dix D) for the ratio (aij/coj)^ equal to 4 and was plotted on 
the same figure. The values given by equation (C2) are 
seen to be in excellent agreement with the theoretically 
exact values. 

The effect of the coupling between bending and torsion 
is to reduce the first natural frequency below that which 
would exist if there were no coupling. This effect is shown 
in figure 4 , wherein the value of 0 is always less than 1. 
This decrease in the first natural frequency due to coupling 
is, however, relatively unimportant in the practical range of 
(co,/co 6)^>4 and 6<C0.75. 

SUMMARY OF RESULTS 

A method based on the use of Station Functions is pre- 
sented for calculating uncoupled and coupled bending-torsion 
modes and frequencies of arbitrary continuous cantilever 
beams. The results of calculations made by this method 


indicated that by the use of Station Functions derived herein, 
n modes and frequencies can be obtained with sufficient ac- 
curacy by using just n stations along the beam if the beam is 
uniform. For a tapered beam, more stations may be re- 
quired, depending on the amount of taper. The amount of 
computational labor is markedly less than for other methods. 
The use of Station Numbers tabulated herein further re- 
duces the amount of calculation necessary. The effect of 
coupling between bending and torsion is shown to reduce the 
first natural frequency to a value below that which it would 
have if there were no coupling. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 18, 1949. 


APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 

dij 

coefficient in equation for Station Function 
in torsion 

B 

bending stiffness of beam, function of z 

Bo 

bending stiffness at root section of beam 


ratio of average bending stiffness of k'’*' 
interval to bending stiffness of root 
section 

bij 

coefficient in equation for Station Function 
in bending 

C 

torsional stiffness of beam, function of z 

Co 

torsional stiffness of root section of beam 

C, 

ratio of average torsional stiffness of P’'' 
interval to torsional stiffness at root 
section 

Cl, C2, C3 

constants defined in appendix B 

/i(2) 

Station Function in torsion for j*'" station 
(defined in text) 

yj{z) 

Station Function in bending for Y" station 
(defined in text) 

I 

mass moment of inertia per unit length of 
beam about elastic axis, function of z, 
except where otherwise defined 

/o 

mass moment of inertia per unit length of 
beam about elastic axis at root section 

h 

ratio of average mass moment of inertia 
per unit length of interval to mass 

moment of inertia per unit length at root 
section 

i,j,k,n 

station indices 

j,k,r 

summation indices 

1 

length of beam 

Mjk, Njlc, P jlC! 

Station Numbers (defined in text); function 

t?,„ N', 

P'm Q'li. 

of indices j, k, and n 

m 

mass per unit length of beam, function of z 

mo 

mass per unit length of beam at root section 

m* 

ratio of average mass per unit length of 
interval to mass per unit length at 
root section 

n 

number of stations along beam 


2o(z) 

2<(3) 




S 

So 

S, 


X 

y 

Vi 

z 

®i7i Ytjj 

r 

T 

s 


6 

0 


0t 

X 




O) 

“i. 


W/ 


bending loading function on beam 
torsional loading function on beam 
absolute magnitude of projection of distance 
from clastic axis to center of gravit}’’ on 
perpendicular to bending direction 
radius of gyration about clastic axis at 
root section 

absolute magnitude of projection of distance 
from elastic axis to center of gravity on 
perpendicular to bending direction for 
root section 

static mass unbalance, function of z, mr 
static mass unbalance at root section, MoTo 
ratio of average of static mass unbalance at 
k^'^‘ section to static mass unbalance 
at root section 

distance fi’om root of beam, except where 
otherwise defined 
bending deflection, function of z 
bending deflection at f*'" station 
dimensionless distance along beam, x/d 
elements of dynamical matrix defined in text 

Co mo 

uncoupled frequency ratio, (w, l^bY 
length of interval along beam between 
two stations 

coupling coefficient, (ro/r^o)^ 
torsional deflection, function of 2 
torsional deflection at i*” station 
root of frequency equation or characteristic 
root of dynamical matrix 
frequency ratio, (u/cos)' 
frequency of vibration 
frequency of uncoupled fundamental bend- 
ing mode 

frequency of uncoupled fundamental tor- 
sional mode 

second derivative of deflection with respect 
to time 
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APPENDIX B 


STATION FUNCTIONS AND DETEKMINANTAL EQUATIONS 


TORSIONAL VIBRATIONS 

A schematic diagram of a cantilever beam divided into n 
intervals of length 5 is shown in figure 1. The Station 
Functions for the torsional vibrations of such a beam must 
satisfy the following conditions; 


2 = 0 

//0)=0 

(Bl) 


/'^(^^)=o 

(B2) 

z=i 

//f)=l 

(B3) 

zH 

/i(j) = 0 j^i 

(B4) 


where/' ( 2 ) denotes the derivative with respect to 2 . 

Equations (Bl) and (B2) represent the boundaiy condi- 
tions that must be satisfied by a cantilever beam vibrating 
in torsion; equations (B3) and (B4) represent the further 
conditions imposed upon the Station Functions. Those 
conditions will be satisfied by a function of the type 

Jii.z) = a\iZ-\-aoiZ^-\- . . . (B5) 


where the coefficients atj must satisfy the following simul- 
taneous equations obtained from conditions (B2), (B3), 
and (B4): 

0 = aii-|- 2 na 2 i + 3n^a3i-|- . . . -/(n-j- (B2a) 

1 =ia,i + 'i^a 2 i + i^a3i+ • • • m (B3a) 

0 = iaii-b.fa 2 i + /«3i+ • • • + jVi (B4a) 

The coefficients a,, can be obtained by solving equations 
(B2a) to (B4a) and the functions ji{z) determined for each 
station. Equation (B5), however, can also be written in 
the following form: 

\l{z—j)z{z—c,) 

where II represents the product for all values of j except 

jy^i 

j=i. The function in equation (B5a) obviously satisfies 
conditions (Bl), (B3), and (B4) because it has zeros at all 
points specified by equation (B4), it equals 1 at the point 
specified by equation (B3), and it equals zero at the point 
specified by equation (Bl). In order to satisfy condition 
(B2), the constant Ci is determined by substitution of equa- 
tion (B5a) into equation (B2). 

Ci=n for 19^71 


Ci=n 


(' 


1 + S 

jF^n 



for i — n 


Equation (B5) can be obtained from equation (B5a) by 
carrying out the indicated multiplications. The complete 
deflection function is then given by 

■ ■ ■ +fniz)0u 
J=l 

The continuous loading function qt ( 2 ) can now be written 
as 

qt{x) = Iw‘^e{z) = Iui^'^Jj{z)ej (B7) 

J=1 

A continuous loading function, which is a function of the 
deflections at the reference stations, has thus been obtained. 

BENDING VIBRATIONS 

The Station Functions for the bending vibrations of the 
beam shown in figure 1 must satisfy the following conditions; 
at 


2 = 0 

p/0) = 0 

(B8) 

2 = 0 

<;'i(0)=0 

(B9) 

z = n 

ry",(7),) = 0 

(BIO) 

z = n 

g"'i{n) = 0 

(Bll) 

z = i 

iyi(f)=i 

(B12) 

z=j 

0iij) = 0 j^i 

(B13) 


where g' (z), g" (z), and g"' ( 2 ) denote the first, second, and 
third derivatives, respectively, of g ( 2 ) with respect to 2 . 

Equations (B8) to (Bll) represent the boundary condi- 
tions that must be satisfied by a cantilever beam vibrating 
in bending and equations (B12) and (Bl3) represent the 
additional conditions imposed upon the Station Functions. 

These conditions will be satisfied by fimctions of the type 


</i( 2)=&2(Z^+63I2^+ • • • (Bl4) 

where the coefficients must satisfy the following equa- 
tions obtained from conditions (BlO) to (B13): 

0 = ‘2b2t-\-&nh^i-\- . . . -)-(?iH“ 3)(?i -1-2)71/^"'*''^ 6 (n+ 3 )i (BlOa) 

0 = 653(-j-247i.64,-l- . . . -i-(7i<-l-3)(7i-l-2)(7i-i- 1 )t!/"6 („+ 3 )i (Blla) 

1 =7^f>2iT'f^f*3iT • • • (/i + 3 )f (Bl2a) 

0 =i^ 62 i+/ 63 i+ • • • -fi‘”'^®*6(n+3)i jVf (Bl3a) 


The coefficients can therefore be obtained from equations 
(BlOa) to (Bl3a) and the fimctions gi ( 2 ) determined for 

11 
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each station i. Equation (B14) can, however, be written 
in the following form: 

U{z—j)Z%Z^ + C2Z+C3) 

where II represents the product for all values of j except 

j~i. The function in equation (Bl4a) obviously satisfies 
conditions (B 8 ), (B9), (Bl 2 ), and (Bl3), because it has 
zeros at all points specified by conditions (B 8 ), (B9), and 
(Bl 3 ) and equals 1 at the point specified by equation (B 12 ). 
In order to satisfy conditions (BlO) and (Bll), the constants 
C 2 and Cs are determined by substitution of equation (Bl4a) 
into equations (BlO) and (Bll). The general forms for C 2 
and C 3 arc, however, complicated and it is easier to obtain 
the numerical values of these constants for each specific 
case. Equation (B14) can then be obtained from equation 
(Bl4a) by carrying out the indicated multiplications. The 
complete deflection function is then given by 

yiz)=YhgAz)yi (Bis) 

j=> 

The continuous bending loading function g^fz) can now be 
M^ritten as 

g&(2) = mwhj(z) = rruo^Yhy A^)y i (B 1 6) 

COUPLED BENDING-TORSION VIBRATIONS 

The Station Functions for the coupled bending-torsion 
vibrations are the same as previously given for the bending 
vibrations and the torsion vibrations. The loading func- 
tions, however, are given as follows (reference 7): 

qt{z) = Ioi^e{z)ArS or y{z) 

= [/ iAz) dj + Sgj(z) yj] (B 1 7) 

j=i 

and 

qi(z) = Sw'^e{z)-]-mu^y(z) 

= [SJj{z)Sj + mgjiz) yj] (BI 8 ) 

j=i 

determinantal equations and dynamical matrices 


0 ,= S^JJg,(z) ^ dz + 8^j\,{z) ^ dz (B19) 

If O is assumed to have a constant value for each interval, 
these integrals may be written as the sum of integrals over 
each section. Equation (B19) then becomes 

= S j] zgtiz) dz^- ^^(l--'k)qi{z) dz-\- ^ ^ qt{z)d^ 

(B20) 

By substituting the relation 


qt{z) = orI '^Jj{z)dj 
;=i 

and by assuming a constant value for I for each interval and 
changing the summation order, 


8i = o} 


Let 


5^^" i: ! s ^ r r dz-ik- 1) /, f mdz+ 

(o; = l (A-=l(ArL Jk-l J k -1 

(B21) 

1 


S 4 f /;(2)dsl|e; 

r=k+\ Jr-1 J) 

i zfAz)dz=Njic 

J k~l 

j^_Uz)dz^M,, J 


Then 


where 


6i~oi^ 7=T “O' 

(■0 j=i 


(B22) 


(B23) 


i] ^ \ h Nj,- ik-1) I,Mj,+ i) Ir Mjrl (B24) 

A-=l Ct L r=t-H 


If C\=It=l (constant cross section), then 
_L r 


Let 


Then 


Tai,., - (* - 1 ) Mj, + 2 Mjrl (B25) 

A — 1 L r — A"j“l J 

^-0 






(B26) 


\9i = '^aij 6j 
j=l 


(B23a) 


Once the Station Functions and the corresponding loading 
functions have been determined, the deflections at the 
reference stations can be obtained in terms of the loading 
function. A homogeneous equation in the reference-station 
deflections for each station is thereby obtained. The 
determinant of the coefficients of the resultant set of homo- 
geneous equations can be set equal to zero; the, determinantal 
frequency equation is thus derived. The deflections at the 
reference stations are obtained by the well-known equations 
for obtaining influence coefficients. 

Torsion. — The deflection at the station i due to the 
continuous loading g,(z) on the beam is given by 


and the characteristic equation is 

IKd-A/| = 0 (B27) 

where I is the identity matrix. 

Bending. — The deflection at the station i due to the con- 
tinuous loading on the beam will be given by 

Jo' d^xdz (B28) 
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If B is assumed to have a constant value for each interval, 
these integrals may be written as the sum of integrals over 
each interval. Equation (B28) then becomes 

*'" 1 ; [j-»- 1)^+5 <*“!)■] 

Jk-I {k — iyz+^ (A; — 1 )“J qi,{ 2 )dz + 

^Ik + 

<l,(z)dz]^ (B29) 

By substituting the relation 


qt,{2) = oihn'£^(jj(z)yj 
J=i 


(B30) 


and by assuming a constant average value for ?n in each in- 
terval and changing the summation order, 


y i 1 } PuVj 

-Oq j=l 


(B31) 


where 




+ 


( 


1 )^ (2k— 1) . 
3 “ Y 




(B32) 


Ja -1 + i (*-l)^] gAz) dz 

' ik= i zgAz) dz 

J k-1 

'«•= f ihiz) dz 

Jk-\ 


^ (B33) 


N 

M 


For a uniform beam, m,^=Bk='\- and equation (B32) be- 
comes 


(iPAk~QAk+±^ {(i-*+|)ATV + 




A:^— (fc — 1)^ (2A: — 1) 


;] Mv|) 


Let 


X = 


Bn 


(jT&^rno 


(B32a) 

(B34) 


then the characteristic equation becomes 


where I is the identity matrix and is the dynamical ma- 
trix. In expanded form, equation (B35) becomes 


dn- 

X dl 2 

• • din 

d 21 

d 22 — ^ • 

• d 2 n 

Pnl 

dre 2 • • 

• dnn — X 


(B35a) 


where X is a latent root of the matrix [/Sy]. 


Coupled bending-torsion vibrations. — The deflections at 
station i are given as before by equations (B19) and (B28). 
The loading functions g, and are changed as follows; 


g,(2) = to-[/ ( 9 ( 2 ) -t- S' 1/(2)] i 
qAz) = (A^[S 0(2) -I- m 2/(2)]) 


(B36) 


If these two equations are substituted into equations (Bl9) 
and (B28) and the integrations arc performed as previously, 
the following relation is obtained: 


~B~ 

>' 0 Po 


1=1 \ > 0 / 

J =1 \ > 0 / 




(B37) 


where a , 7 and ,3 , 7 are given in equations (B24) and (B32) and 


rgo 


1 hBo 

'o'f'iO 


^ S- C\:m 


Tu=S y-t r SkN' jic—{k — ]) S/cM' jkA- S, 

A-=l La-L = J 

5.7'=S Sk[i Pjk— Qjk]-\- 

i) j (B38) 

where 

Pjk^ Ja -1 (/r- 1 )^] { 2 ) dz 

Qik = [^^“1 (A;— 1)^2 +1 {k—lf~^Jj{z)dz 

the determinantal equation therefore is 


[voll = 0 


where [vij] is tbe dynamical matrix, the elements of which 
are as indicated in equation (B37). The matrix [ 17 ^] is seen to 
be a 2 tiX 2 n matrix. 


l[(3y]-X/l = 0 


(B35) 


APPENDIX C 


QUADRATIC FORMULA FOR FIRST COUPLED MODE 


If only the first vibrational mode is desired, it is possible 
to obtain this mode approximately by coupling together the 
fundamental uncoupled bending mode with the fundamental 
uncoupled torsional mode to obtain a simple quadratic 
equation for the first coupled frequency. This equation is 
valid when the coupling coefficient e is constant along the 
beam. The differential equations obtained by coupling the 
fundamental uncoupled torsional mode with the funda- 
mental uncoupled bending mode are: 

m V H- S e + m coj- V ~ 0 ) 

> (Cl) 

Sy Id 1 / 

where 

m mass per unit length of beam, function of z 
S static mass unbalance, function of z 


I mass moment of inertia about elastic axis, fimction of z 
03„ frequency of uncoupled fundamental bending mode 
CO, frequency of uncoupled fundamental torsional mode 
. . denotes differentiation twice with respect to time 

These equations lead to a quadratic equation in the fre- 
quencj' ratio fi, whose solution for the lowest frequencj^, 
provided e is constant along the beam, is 


where 


fi=- 


COs 


1—7 
2(1 -e) 



4t(1 — tf 


(C2) 


0 frequency ratio, (co/co^)- 
7 uncoupled frequency ratio, (co,/coj)^ 

€ coupling coefficient, (r//'g)^ 

This quadratic has been plotted in figure 4 for values of e 
ranging from 0 to 1 and values of 7 = (co, /coj)^ from 1 to 100. 


APPENDIX D 

EXACT SOLUTION FOR COUPLED BENDING-TORSION VIBRATIONS OF UNIFORM CANTILEVER BEAM 

where 


The differential equations for the equilibrium of an ele- 
ment of a beam vibrating in coupled bending-torsion vibra- 
tions can be put in the following dimensionless form: 


d‘* Fi ml^ 


dx* 

d^Y2 


B 




IV 


B 

IV 

2 V , ,2 V 

j I — ^ u J 2 


where 


X =- 


Now 


Y.^ylr 

Y2=e 

distance from root 
I 


, CiB 


2 _ 

^5 JJ2 


where 


C4=12.36 

C5=2.467 

Equations (Dl) become 


- e X 1 - 
7 


dxV 




Y 2 


(Dl) 


(D2) 


Let 


Then 


dY, 

dx 

dYl. 

dx 

dY, 
dx ■ 

dYn _ 

dx 

dY, 

dx 

dY, 


n = (co/c06)- 
7 = (c^,/co6)^ 

= Y, 

-Y, 

- CjOl } 1 1 2) 

(,Y,+ Y2) 


(D3) 


dx 


Equation (D3) can be written as the single mati’Lx equation 


A 

dx 


"IT 


~ 0 

0 

1 

0 

0 

0" 


"IT 

Y 2 


0 

0 

0 

0 

0 

1 


Y 2 

F 3 


0 

0 

0 

1 

0 

0 


Y, 

Y, 

= 

0 

0 

0 

0 

1 

0 


T4 

Y, 


C 4 n 

C 40 

0 

0 

0 

0 


F5 

Y, 


1 

1 

— Csfl 
7 

0 

0 

0 

0 


Fa 


(D4) 
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or 


dx 


= AY 


(D4a) 


The X matrix of A is 

r — X 0 1 


0 0 


0 1 


where Y and A are the matrices indicated. 

The solution to the matrix equation (D4) is given by 

F=e^"Fo (D5) 


where Fo is a column of arbitrary constants. 
From the boundary conditions 

at x = 0 Fi=F 2 =F 3=0 

x=l F4=F3=Fe=0 


F„=F(0) = 


0 ■ 
0 
0 

F4(0) 

, 1^5(0) 

LF(o). 


If then is an element of the matrizant the boundary 
conditions give 


0 

— X 

0 

0 

0 

1 

0 

0 

— X 

1 

0 

0 

0 

0 

0 

— X 

1 

0 


C4I2 

0 

0 

— X 

0 

C5S2 

y 

Cfll2 

y 

0 

0 

0 

— X 


The characteristic equation A(X)=0 is 

X«+^ X^-C4fiX2-(l -e)c4C5 Y=0 (D7) 

Equation (D7) is a cubic equation in X^. Let the roots be 

Xi, Xi, X2, — X2, X3, — X3 

Then by the confluent form of Sylvester’s theorem, 


fi44 




1245 


= 0 


(D6) 


I 1254 Qgj 1250 I 

Equation (D6) is the frequency equation. It has an infinite 
number of roots for co. 

In order to determine the elements 12;^, must be evalu- 
ated. Use will be made of Sylvester’s theorem (reference 
13). 


1 r g^^(x) -| 

l(ai-i) dX“<-i I n(x-x*)“*^ 

Jx=x 


(D8) 


where F{\) is the adjoint matrix, r is the number of distinct 
roots, and at is tbe multiplicity of tbe i'" root. 

If the roots are all distinct, this relation becomes 


2X4n(Xi-x,)(Xi+x,) 

jVi 

where the adjoint matrix F(X) is given by 


(D9) 


X=-f — x^ 
y 

C 4 fiX 

x 4 +£^ x= 

7 

X'+^X 

7 

X'^+Cj- 

7 

C 4 I 2 

-eC 3 ^X^ 

T 


-€^X^ 

7 

C 5 I 2 .. 
— € X 

y 

Csl 2 

— 6 

y 

X^ — C 4 O 

02 

C40X^ + (1 — e)c 4 Cs — 

y 

C412X2 

x®-t-^ X" 

7 

xi+£^ x^ 
7 

X3+^X 

7 

C 4 I 2 X 

02 

C40X^-t-(l — €)C4C5 — X 

y 

C4l2X^ 

C4nx^-l-(i «) 

7 

x=-t-^ x^ 
7 

X^-b^ X^ 
7 

C4$2X2 

C 40 X‘- 1 -( 1 — e)c 4 C 0 — X" 
7 

C4l2X‘ 

C4i2X^+(l t) X 

y 

C4l2X2-l-(l-e) 

y 

X«-b— X® 
7 

C 4 OX® 

eCsll 4 
7 

7 7 

-e^X^ 

7 

- 6 ^x 2 

7 

-e^X 

7 

X'— C412X 


(DIO) 


16 


RBPOBT 1005— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


From equations (D9) and (DlO), the elements Qa are seen 
to be given by 




\ 4 I X 2 

3 


Mi 


cosh Xj 


^^45 ^ ■ 


3 


Mi 


sinh Xj 


3 C40Xi" . , , 

Mi 


o\ 2 I ^4^6^^ /-I \ 

:fiXi^H (1— e) 


H54 


f ,r x;n(Xi2-x/) 

iT^i 


sinh Xi 


O^s — ^^4 


y (Dll) 


“‘•=-Sn(x7=rS“*'' 

j^i 


eCs \i 

^ sinh X i 


0^4 — ■ 


Iti n(Xi^-x/) 

Mi 

-€C--X 

XiH (xAx^^ 

jVi 

C 4 S 2 1 , 

^-=-§ 

Mi 


The value of the determinant in equation (D6) must be plot- 
ted against the frequency; the value of the frequency for 
which this determinant becomes zero is thereby obtained. 
This procedure involves first solving the cubic equation (D7) 


for each assumed value of frequency parameter and then 
calculating the elements of the determinant from equations 
(Dll). The process is evidently long and laborious. 
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TABLE I— STATION NUMBERS 
n=l 


TABLE II— STATION NUMBERS 
n = 2 



k 

j 

1 

2 



11 

5 

M 

1 





12 

12 

N 


8 

8 



15 

15 

P 


0. 183333 

0. 025000 

Q 


. 046032 

. 029365 

M' 


. 536364 

. 627273 

N' 


. 367100 

. 851948 

P' 


. 137933 

. 057955 

Q' 


. 036616 

. 069733 



13 

29 


2 





48 

48 



31 

239 

N 


^0 

240 

P 


-0. 037500 

0. 143750 

0 


008135 

. 181448 

M' 


— . 060795 

. 448674 

N' 


-.034875 

. 758685 

P' 


-.011252 

. 118462 

Q' 


-. 002614 

. 150415 
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TABLE III— STATION NUMBERS 


n = Z 



J 

1 

2 

3 

M 

1 

0. 950000 

0. 450000 

-0. 050000 

N 


. 645833 

. 587500 

-.120833 

P 


. 186310 

. 032143 

005357 

Sj 


. 046577 

. 038244 

-.011756 



. 596268 

. 533205 

097426 

N' 


.390646 

. 708205 

239994 

P' 


. 148013 

. 042560 

-. 012798 

Q' 


. 038884 

. 050843 

-. 028318 

M 

2 

-0. 525000 

0. 725000 

0. 475000 

N 


24I6G7 

1. 175000 

1. 091667 

P 


068452 

. 160714 

. 031548 



— . 014583 

. 202083 

. 068750 

M' 


149356 

. 602896 

.625418 

N' 


083406 

. 994860 

1. 475153 

P' 


026378 

. 143948 

. 057937 

Q' 


006034 

. 181698 

. 127659 

M 

3 

0. 235185 

-0. 153704 

0. 568519 

N 


. 106019 

231944 

1.513426 



. 029563 

— . 023677 

. 139749 



. 006222 

028963 

. 316408 

M' 


. 040630 

072928 

. 445812 

N' 


. 022325 

111744 

1. 200133 



. 006972 

-.012081 

. 118007 

Q' 


.001579 

014830 

. 267865 


TABLE V— STATION NUMBERS 


n = 5 




1 

2 

3 

4 

5 

M 

1 

1.097991 

0. 408755 

-0. 040898 

0. 019866 

-0.013120 

N 


. (>08222 

. 527493 

-. 100112 

. 069159 

—.058445 

P 


. 202887 

. 02(5908 

-.005074 

. 002776 

— .001627 

Q 


. 049943 

.031910 

-.011210 

. 008927 

— . (H16836 

M' 


. 649902 

.492141 

-.070298 

. 034939 

-.024007 

N' 


. 427616 

. 647530 

-. 172488 

. 121619 

— . 107639 

P' 


. 156411 

. 03(5908 

— . 008903 

. 004824 

— . 00.3:375 

Q' 


. 040729 

. 043977 

-.019678 

. 015509 

-.014228 

M 

2 

-0. 839550 

0. 799339 

0. 49.3783 

-0. 089550 

0. 049339 

N 


-.373049 

1.282044 

1. 145470 

-.310549 

. 219544 

P 


-. 103119 

. 169599 

. 037952 

-.011949 

. 006007 

0 


-.021583 

, 212792 

. 083245 

-.038398 

. 02524:3 

M' 


-. 255330 

. 699256 

. 523828 

-.09972:3 

. 0581:34 

N' 


-. 139170 

1. 138472 

1.219101 

-.345551 

. 260447 

P' 


-.043239 

.157833 

.041704 

-.013166 

. 008078 

Q' 


-.009758 

. 198610 

.091532 

-.042299 

. 0:34055 

M 

3 

0. 762798 

-0.313591 

0. 651687 

0. 575298 

—0. 126U91 

iV" 


.329315 

-. 465823 

1.718204 

1.923065 

— . 559573 

P 


. 089079 

-.044602 

. 160488 

. 049347 

—.014691 

0 


. 018334 

-.054326 

. 340126 

. 157843 

-. 061693 

iWf' 


. 197103 

-. 228783 

. 633812 

. 573549 

— . 137821 

JSf' 


. 105126 

-.344915 

1.674943 

1.916360 

-.6162.34 

P' 


. 032115 

-.034985 

. 148520 

. 048803 

-.018600 

Q' 


. 007151 

-.042759 

. 335798 

. 156076 

-. 078385 

M 

4 

—0. 548214 

0. 187897 

-0. 159325 

0. 576786 

0.562897 

N 


-.233780 

. 276637 

-. 400446 

2. 109970 

2. 432887 

P 


-. 062665 

. 025479 

-.024868 

. 140460 

.042295 

0 


-.012809 

. 030950 

-.065302 

. 458397 

. 177006 

Af' 

"KT! 


— . 117990 

. 117132 

-. 136452 

. 557359 

.664560 

JS 


— . 062435 

. 175188 

-.344108 

2. 041363 

2. 901646 

P' 


-.018958 

.017186 

-.021868 

. 137234 

. 063512 

Q' 


-.004201 

. 020954 

-.048664 

. 447980 

. 267043 

M 

5 

0. 214238 

-0. 069026 

0. 050904 

-0. 080137 

0. 525349 

N 


. 090928 

-.101352 

. 127410 

-. 283759 

2. 458336 

P 


. 024289 

009225 

. 007667 

-.013645 

. 134478 

0 


. 004952 

-.011196 

.017031 

-.044065 

. 573757 

Af' 


. 033722 

-.031711 

. 032307 

-. 056459 

. 432107 

N' 


. 017789 

-.047311 

.081140 

- 200039 

2. 030260 

P' 


. 005389 

-. 004593 

. 005002 

-.009675 

. 116059 

Q’ 

1 . U01192 1 

— . 005596 

.011120 

-.031249 

.495663 


TABLE IV— STATION NUMBERS 


n = 4 




1 

2 

3 

4 

M 

1 

1. 022222 

0. 429630 

-0. 051852 

0. 022222 

N 


. 576455 

. 557937 

127249 

. 076455 

P 


. 194478 

. 029597 

-.006581 

. 002612 

Q 


. 048240 

. 035167 

-.014547 

. 008359 

M' 


. 623188 

.511882 

-.082891 

. 042276 

N' 


. 413738 

. 676680 

203719 

. 146954 

P' 


. 152256 

. 039616 

-.010661 

.005795 

O' 


. 039818 

. 047267 

-. 023551 

.018630 

M 

2 

-0. 647917 

0. 747917 

0. 618750 

-0. 085417 

N 


-.292857 

1. 207143 

1. 207143 

292857 

P 


-.081920 

. 163021 

. 041295 

-.009598 

O 


-. 017295 

. 204828 

.090642 

-.030688 

M' 


-.211987 

. 667412 

. 544025 

-. 112648 

N' 


-.116602 

1.091462 

1. 269193 

-.390685 

P' 


-.036502 

. 153469 

.044508 

-.015000 

O' 


—.008281 

. 193310 

. 097745 

-.048203 

M 

3 

0. 522222 

-0. 255556 

0. 633333 

0. 522222 

N 


. 229365 

-.381746 

1.673810 

1.729365 



. 062798 

-.037401 

.148512 

. 037202 

0 


. 013040 

-.045624 

. 335791 

.118397 

M' 


. 122052 

-.166738 

. 582158 

. 643846 

N' 


. 065879 

252823 

1. 545802 

2. 164827 



. 020304 

-.026235 

. 140822 

. 060554 

0' 


. 004551 

-.032114 

. 318707 

. 194016 

M 

4 

-0. 221701 

0. 094850 

-0. 105961 

0. 543924 

N 


096544 

. 140724 

-.267841 

1.997206 

P 


-.026267 

.013391 

-.017322 

. 136803 

0 


-. 005428 

.016301 

-.038574 

. 446729 

M' 


-. 035456 

. 042628 

064723 

. 438962 



-.019023 

. 064205 

164169 

1. 622066 



-.005836 

. 006481 

-.010869 

. 117037 

V 1 1 

-.001303 

.007917 

-.024222 

. 382752 


TABLE VI— STATION NUMBERS 

n=& 



A- 

j 

1 

2 


4 

5 

6 

M 

1 

1. 172073 

0. 391101 

-0. 032:171 

0. 01332:1 

-0.010149 

0. 008879 

N 


. 638800 

. 5018.56 

-. 078978 

. 046300 

-. 045644 

. 048522 



. 210893 

. 024685 

003894 

. 001823 

— . 001498 

.00114.3 

Q 


. 051551 

. 029221 

008595 

. 005860 

006322 

. 005949 

M' 


. 676394 

.474177 

-.059129 

. 026582 

-.018685 

. 015649 



. 441269 

. 621067 

-. 1447.59 

. 092370 

-. 083901 

. 085903 



. 160476 

. 034473 

007337 

. 003631 

002691 

. 002241 

Q' 


. 041616 

. 041021 

-.016206 

. 011674 

011350 

. 011694 

M 

2 

-1.066598 

0. 853106 

0. 468124 

-0. 070505 

0. 044513 

-0. 035782 



466718 

1. .3()()105 

1. 081893 

— . 24 1062 

. 19994H 

— . 195451 



-. 127634 

. 1 763.58 

. 034411 

009203 

. 006452 

— . 004561 

0 


026.505 

. 220969 

. 075401 

029565 

.027210 

023740 



.30.3948 

. 731991 

. .503742 

-.085145 

. 049793 

— . 038661 



-. 164215 

1. 186681 

1. 169252 

2947.36 

. 223322 

— . 212149 



-. 050692 

. 162263 

. 038896 

011103 

. 007043 

005503 

Q' 


-.011384 

. 203987 

. 085309 

-. 035665 

. 029701 

-. 028708 

M 

3 

1. 150584 

-0. 404200 

0. 693794 

0. .546418 

-0. 133360 

0. 089627 



. 489124 

597296 

1. 822457 

1. 822457 

-. 597296 

.489124 

p 


. 130870 

055956 

. 1 56259 

. 046293 

-. 018484 

.011225 

Q 


.026710 

-. 068060 

. .352909 

. 144808 

-.077925 

. 058410 



.267118 

-. 275585 

. 662165 

. .553477 

126314 

. 083246 



. 141177 

-. 413819 

1. 745286 

1.846431 

564890 

. 456509 



. 042839 

-.041308 

. 152473 

. 045980 

-.017100 

.011713 

Q' 


. 009400 

-. 050434 

. 344556 

. 146998 

072063 

. 061097 

M 

4 

-0. 930965 

0. 273028 

-0. 194854 

0. 502473 

0. 624591 

-0. 171416 



-. 390902 

. 399897 

488124 

2. 163786 

2. 720209 

-.933437 



103C>:}5 

. 036020 

029594 

. 142229 

. 056698 

020561 

0 


021011 

. 043690 

065761 

. 464052 

. 238215 

— . 106938 



210538 

. 182630 

— . 180822 

. 598464 

. 604424 

— . 166643 



110263 

. 271857 

-. 454622 

2. 185427 

2. 628772 

-.912352 



033225 

.026154 

-.028221 

. 143453 

. 053359 

022768 

O' 


007320 

. 031846 

062753 

. 468015 

. 224100 

-. 118729 

M 

5 

0. 581796 

-0. 156399 

0. 092907 

-0. 111954 

0. 537351 

0. 599157 



. 242612 

228221 

.231501 

394888 

2. 509279 

3. 194001 



. 063996 

-.020218 

. 013474 

-.018261 

. 134578 

. 046991 

Q 


. 012925 

-.024496 

. 029896 

-. 058920 

. 574036 


M' 


. 120308 

-.097119 

. 081668 

-.110987 

. 535339 

. 685021 



. 062735 

-. 144101 

. 204039 

-.391731 

2. 499695 

3. 678582 



. 018841 

-. 013674 

. 012212 

-. 018226 

. 133988 

. 066470 

O' 


. 004140 

-.016634 

.027119 

-. 058815 

. 571521 

. 345986 

M 

6 

-0. 209220 

0. 054246 

-0. 030239 

0. 031561 

-0. 064035 

0. 610543 



086982 

. 079042 

-. 075223 

. 110987 

-. 291427 

2. 902746 



— . 022893 

. 009657 

— . 004323 

. 004969 

-.011243 


0 


—.004616 

. 008425 

-. 009587 

. 016022 

-.047575 

. 698254 



033141 

. 025961 

-. 020586 

. 024431 

049677 

.426817 



-. 017248 

. 038471 

—.051417 

. 085977 

-. 225999 

2. 427820 



-.005173 

. 003631 

— . 003042 

. 003877 

— . 008676 


0' 


—.001135 

. 004415 

-. 006753 

. 012502 

036708 

. 606979 
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TABLE VII— STATION NUMBERS 

n— 7 





2 

3 

4 

5 

6 

7 

M 

1 

1.243487 

0.376396 

-0.026266 

0. 009112 

-0. 005896 

0. 006025 

-0. 006513 

N 


. 667840 

. 430602 

-.063889 

.031590 

026481 

. 033195 

-.042160 

P 


.218415 

. 022882 

003069 

.001211 

-. 000853 

. 000925 

-.000363 

0 


. 053049 

. 027042 

-. 006769 

.003890 

003599 

.004829 

-.005357 

M' 


. 702228 

. 458303 

-.050122 

. 019989 

-.012820 

.011370 

-.011099 

A' 


. 454474 

. 597757 

— . 122429 

. 069352 

057533 

. 062529 

-.072080 

P' 


. 1643S2 

. 032358 

-.006088 

. 002679 

001831 

. 001689 

-.001614 

O' 


. 042465 

. 038456 

-- 013441 

.003611 

-.007724 

.008315 

-.010037 

M 

2 

-1.321299 

0. 906437 

0. 446479 

—0. 065674 

0. 029812 

-0. 027896 

0. 028701 

N 


570270 

1. 435730 

1.028397 

-. 192270 

. 133730 

153603 

. 185730 

P 


154452 

. 182773 

.a31489 

007052 

. 004233 

004239 

. 003780 

0 


-.031847 

. 228718 

. 068932 

-.022639 

. 017853 

-.022136 

. 023463 

M' 


-.357636 

. 764868 

. 485193 

— . 071706 

. 038132 

-.031029 

. 028988 

N' 


191650 

1.234950 

1. 123273 

247832 

. 170920 

170555 

. 188216 

P' 


-.058807 

. 166641 

. 036328 

009168 

. 005346 

-.004565 

. 004199 

O' 


-.013146 

. 209297 

. 079621 

-. 029438 

. 022542 

-. 023823 

.026110 

M 

3 

1.672922 

-0. 511106 

0. 737737 

0. 516672 

-0. 104S56 

0. 081487 

-0. 077078 

N 


. 701415 

-.751768 

1.931045 

1. 718603 

-.468955 

. 448232 

-. 498585 

P 


. 185835 

-. 069049 

. 162183 

. 040990 

-. 014225 

.012162 

-.010056 

0 


. 037666 

— . 033877 

.366024 

. 130958 

-.059956 

. 063490 

-. 062408 

M' 


.355047 

-.329159 

. 692136 

.532099 

108454 

.073511 

-.063669 

N' 


. 186093 

492459 

1. 819566 

1.771836 

484634 

. 403671 

413268 

P' 


. 056122 

-. 048434 

.156613 

. 042913 

-.014528 

. 010624 

-.009166 

0! 


.012374 

-.059075 

. 353730 

. 137128 

-.061218 

. 055439 

056987 

M 

4 

-1.605312 

0. 409927 

-0. 250374 

0. 629063 

0. 592219 

-0. 182666 

0. 144688 

N 


-.664780 

. 597643 

— . 625274 

2. 291470 

2. 574726 

-1.002357 

. 935220 

P 


-. 174510 

. 052757 

-.037056 

.147488 

.051989 

-.026098 

. 018551 

0 


-.035122 

.063907 

-. 0S2279 

. 480975 

. 218352 

136173 

.115110 

M' 


-.317567 

. 247432 

-. 216700 

. 623577 

. 584350 

-.157315 

. 115706 

N' 


-. 164912 

. 366953 

— . 543417 

2. 273022 

2. 538692 

-.861833 

. 750618 

P' 


049381 

. 034761 

-. 033166 

. 147041 

. 050504 

-.021773 

. 016464 

O' 


-. 010326 

. 042284 

-.073703 

. 479567 

. 212061 

-.113560 

. 102355 

M 

5 

1. 129029 

-0.264373 

0. 134585 

-0. 136596 

0. 551252 

0. 668960 

-0. 220971 

N 


. 464325 

-.384008 

. 334325 

480675 

2. 570992 

3. 589325 

-1.425675 

P 


. 121269 

-. 033334 

.019010 

-.021700 

. 136202 

. 063551 

-.027154 

0 


. 024312 

-.0403:« 

.042147 

-.069981 

. 580355 

. 330701 

-. 168426 

M' 


. 234133 

— . 168109 

.122950 

— . 143020 

. 568533 

. 634247 

197770 

N' 


. 120971 

-. 248403 

. 306705 

503685 

2. 649577 

3.397336 

-1.281188 

P' 


. 036034 

-.023166 

. 017983 

-.022914 

. 139031 

.057851 

-.027322 

O' 


. 007887 

-. 028149 

.a39907 

-.073907 

. 593022 

. 300912 

-. 169318 

M 

6 

-0.617807 

0. 137401 

-0. 064103 

0. 054068 

-0. 084474 

0. 507772 

0. 632193 

N 


252961 

. 199169 

-. 158887 

. 189539 

-. 383331 

2. 883613 

4. 007039 

P 


-. 065852 

. 017129 

— . 008878 

.003211 

-. 014243 

. 130010 

.051386 

0 


-.013170 

.020712 

019672 

. 026463 

-.060230 

. 684756 

. 318021 

Af' 


-. 125078 

. 086131 

-. 058376 

- 057414 

-. 092005 

. 516450 

. 704705 

N' 


064447 

. 127049 

-. 145336 

.201491 

-.417397 

2. 931074 

4. 491463 

P' 


-.019183 

.011759 

-.003393 

. 008826 

-.015450 

. 131144 

. 069349 

O' 


-. 004186 

- 014281 

— . 018616 

. 028441 

065330 

. 690659 

. 430362 

M 


0. 205449 

-0, 044613 

0. 020059 

-0. 015852 

0. 021382 

-0. 053078 

0. 498306 

N 


. 083943 

-.064603 

.049675 

-. 055505 

. 096798 

-.295079 

3. 334065 

P 


.021819 

-. 005533 

. 002756 

-. 002373 

. 003481 

-.009551 

. 130932 

0 


. 004358 

-.006689 

. 006106 

007644 

.014713 

-. 0499S1 

. 820708 

M' 


. 033023 

-. 022307 

.014642 

-.013579 

. 018939 

-.044219 

. 420133 

N' 


. 016993 

-.032878 

. 036425 

-. 047602 

. 085721 

245636 

2. 816717 

P' 


. 005053 

-. 003a33 

.002091 

002060 

. 003074 

-.007854 

.114318 

O' 


.001102 

-. 003682 

. 004637 

-. 006638 

.012992 

-.041094 

. 716958 
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TABLE VIII— STATION NUMBERS 

n=8 




1 

2 

3 

4 

5 

6 

7 

8 

M 

1 

1.312192 

0. 364019 

-0.021829 

0.006490 

-0. 003545 

0. 003081 

-0. 003931 

0.005039 

N 


.695399 

. 462793 

-.052953 

.022453 

-.015893 

.016958 

-.025621 

. 037691 

P 


.225483 

. 021401 

-.002485 

.000839 

-.000499 

.0Ck)464 

-.000622 

.000684 

9., 


.054447 

.025255 

-.005476 

.002694 

-.002103 

.002422 

-.003869 

.004929 

M' 


. 727233 

.444367 

-.043004 

.015242 

008698 

.007082 

007522 

.008334 

N' 


. 467152 

. 577362 

-. 104816 

.052799 

-.038988 

.038933 

-. 048943 

.062492 

P' 


.168111 

. 030533 

-.006119 

.002003 

-.001220 

.001045 

-. 001145 

.001228 

Q' 


. 043271 

. 036245 

-.011294 

.006433 

-.005146 

.005456 

-.007123 

.008863 

M 

2 

-1.600390 

0. 955663 

0.428499 

-0. 045077 

0. 020351 

-0. 016189 

0. 019610 

-0. 024337 

N 


-.682210 

1. 507998 

.984086 

-. 155335 

.091123 

-.089039 

.127790 

-. 182002 

P 


-.183160 

. 188781 

.029123 

-.005550 

.002808 

-.002412 

.003082 

-.003289 

Q 


-.037524 

.235968 

.063696 

-.017807 

.011838 

-.012593 

.019183 

023705 

M' 


-.415706 

. 797175 

. 468735 

-.060778 

. 028715 

-.021414 

. 021627 

-.023284 

iV' 


-.221089 

1.282229 

1. 082551 

-.209723 

. 128559 

-. 117660 

.140686 

174560 

P' 


-. 067464 

. 170868 

. 034087 

—.007611 

.003953 

-.003129 

. 003272 

-.003421 

Q' 


—.015018 

. 214419 

. 074659 

— . 024428 

.016664 

-.016330 

.020358 

-. 024694 

M 

3 

2. 337500 

-0. 630523 

0. 780389 

0. 491025 

-0. 082648 

0. 054514 

-0. 060500 

0. 071477 

N 


. 967933 

-.923581 

2. 03G155 

1. 631308 

-. 368956 

.299530 

-.394067 

. 534419 

P 


. 254178 

-.083331 

. 167794 

. 037442 

-.010886 

. 007975 

-.009418 

.009602 

Q 


.051181 

101109 

. 378439 

.119546 

-.045861 

. 041626 

-.058605 

.069206 

M' 


. 462586 

-.388839 

.722423 

.512112 

-.091426 

.056607 

-. 052469 

.054058 

N' 


. 240602 

579786 

1. 894495 

1. 702145 

-, 408100 

.310711 

-.341167 

. 405216 

P' 


. 072148 

-.056240 

. 160734 

.040069 

-.012016 

. 008120 

-.007864 

.007911 

Q' 


. 015838 

-.068533 

.362854 

. 127977 

-.050617 

. 042369 

-.048918 

. 057107 

M 

4 

-2. 030070 

0. 593584 

-0.315718 

0. 667195 

0. 658819 

-0. 143453 

0. 132430 

-0. 143916 

N 


-1.075644 

. 861856 

—.786292 

2. 424356 

2. 424356 

-.780292 

. 80 1850 

-1.075644 

P 


-.279850 

. 074709 

-.045646 

. 152881 

. 040981 

-.020073 

. 020253 

-.019152 

Q 


-.055949 

. 090392 

-. 101283 

. 498324 

. 197214 

-.104712 

.126011 

-. 138032 

M' 


-.466670 

. 329868 

-. 258293 

. 650770 

.661588 

135711 

. 104892 

-.100193 

AT' 


-.240468 

. 487559 

-.646305 

2. 367821 

2- 436346 

-.743106 

. 081473 

-.750848 

P' 


-.071591 

. 045534 

-.038818 

.150902 

.047159 

-.018614 

. 015449 

-.014572 

Q' 


-.015627 

. 055337 

-.086223 

. 491981 

. 197950 

-. 097075 

. 096084 

-.105178 

M 

n 

2. 192995 

-0. 454018 

0. 201514 

-0. 175130 

0.584107 

0. 633389 

-0. 237005 

0. 216982 

N 


. 890699 

-. 656778 

. 499217 

-.614926 

2. 718847 

3. 393592 

-1.539301 

1.613222 

P 


. 230544 

-.05597.5 

. 027793 

-.027134 

. 141097 

. 058226 

-.034728 

. 028253 

Q 


.045912 

067648 

.061577 

-.087460 

.601511 

. 302906 

-.215987 

. 203600 

M' 


. .389428 

-.253524 

. 165704 

— . 170613 

.591124 

.614158 

191506 

. 155413 

N' 


. 199636 

-.373346 

.412443 

599857 

2. 750891 

3. 287044 

-1.241724 

1. 164103 

F' 


. 059201 

-.034316 

. 023783 

-.026822 

.142374 

. 054969 

-.027027 

. 022340 

<?' 


. 012882 

-.041657 

. 052749 

-.086477 

. 006915 

. 285874 

168026 

. 161236 

M 

6 

-1.350286 

0. 265559 

-0. 108064 

0. 078152 

-0. 102878 

0. 520374 

0. 709714 

-0. 274441 

N 


-.546005 

. 383370 

-.207117 

. 273370 

-. 466005 

2. 952258 

4. 523995 

-2. 046630 

P 


-. 140865 

. 032367 

-.014611 

.011568 

016907 

. 131525 

. 070024 

-.034431 

Q 


-.027983 

. 039093 

0;i2351 

. 037251 

-.071468 

. 692634 

. 434540 

-.248035 

M' 


-. 263362 

. 163432 

-.098083 

.084195 

-.115917 

. 543778 

. 662638 

-.230699 

N' 


-. 134575 

.240192 

-.243600 

.294907 

525024 

3. 081687 

4. 216400 

-1.725908 

P' 


-.039810 

.021885 

-.01.3810 

. 012656 

-.019020 

. 135367 

. 062213 

-.032195 

Q' 


-.008640 

. 026552 

-.030611 

. 040765 

-. 080395 

.712714 

. 385904 

-.232314 

M 

7 

0. 654484 

-0. 124395 

0. 048081 

-0. 031856 

0. 034801 

-0.066829 

0. 484484 

0. 662747 

N 


. 203837 

-. 179348 

. 118705 

111252 

.157170 

-.370670 

3. 239551 

4. 867795 

P 


. 067914 

-.015052 

. 006430 

004626 

.005472 

-.011570 

.126311 

. 055545 

Q 


.013468 

-.018172 

.014232 

-.014889 

.023114 

-.060517 

. 791745 

. 399381 

M' 


. 130945 

-.079057 

.045388 

~. 036054 

.041862 

-.077799 

.500211 

. 723527 

N' 


. 066783 

-. 116059 

. 112600 

-. 126105 

. 189085 

-.431248 

3.341811 

5. 337321 

P' 


. 019727 

-.010524 

. 006331 

-.005328 

.006593 

-.013324 

. 128649 

. 072134 

Q' 


. 004280 

-.012763 

.014029 

-.017155 

. 027846 

-.069680 

. 806281 

. 519815 

M 

S 

-0.202414 

0. 037821 

-0.014268 

0. 009104 

-0. 009307 

0.015388 

-0. 045167 

0. 487926 

N 


-.081470 

. 054494 

035205 

.031774 

-.041992 

.085181 

-.296574 

3. 754721 

P 


-.020947 

. 004560 

— . 001898 

. 001312 

-.001443 

. 002573 

-.008294 

. 129520 

Q 


-.004150 

.005504 

004201 

. 004221 

-.006093 

.013452 

-.051711 

. 941457 

M' 


-.033107 

. 019722 

-.011093 

.008527 

-.009326 

.015058 

-.039785 

.414996 

N' 


-.016869 

. 028937 

027507 

.029809 

-.042089 

. 083310 

- 260993 

3. 198326 

P' 


—.004979 

.002018 

-.001541 

.001252 

-.001450 

. 002493 

-.007173 

. 113558 

Q' 


-.001080 

.003174 

-.003414 

. 004029 

-.006122 

. 013034 

-.044710 

. 825790 
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TABLE IX— COMPARISON OF RESULTS 


Numberl 

of 

stations 


Torsion 


Bending 


Coupled 


Wj 


/of 

Co 


'2-y/- 


Bo 





Station-I 

'unction raetho< 

i 




1. 549 



3. 493 



3.46 



2 

1. 571 

4. 520 


3. 516 

21. 71 


3. 48 

19. 7 


3 

L571 

4. 689 

7.502 

3. 516 

22.04 

60. 20 

3. 48 

20.6 

48.2 

Weighted influenoe coeflicicnts 



5.39 


3. 56 

15. 63 






1.571 

4. 73 


3. 52 

22.80 














Exact theoretical value 


1.571 

4.712 

7. 854 

3.516 

22.04 

61.70 

3. 49 

20.6 

49. 1 


TABLE X— STATIONS REQUIRED FOR SATISFACTORY ACCURACY 


Method 

Torsion 

Bending 

V Co 

V Co 

V Co 

V Bo 

V Bo 

V Bo 

Station Functions 

1 

3 

4 

1 

2 

3 


2 

4 


3 

6 










U. S. government printing OFFICE: 1951 




V ^ / 


■ ^ :j 




V •. 




y" 



y Positive directions of axes apd angles forces an<i moments) -aite shown by- arrows 




-■ J . 1 


' y Axis 

- 

.^Force- 
(parallel 
to axis) 
symbol 

■i' ' ■■ ^ -1 - ' 

Moment- about axis r 

Angle 

y:- 

' Velocities . ^ 

-S 

Designation 

Sym-- 
bol y 

N"" 

^Designation 

7<- ) - 

1 • - 

Sym- 

bol 

' Positive ; 
-'-ilireotibh^ 

iDesigna^- 
Y tion ^ 

Sym- 
jbol ' 

'Linear 
(compo- 
nent along 
axis)- 

Angular 

LongitudfnaL r 
. Lateral...'._-i.'... 

X 

y- 

X ■ 

- ~ r 

JRoUi^„..ll 
Pitching., .. 

t-L '■ 

M 

Z^-r^X-^ 

BoU-... 

Pitnh 

« ' 
\0 

■■ U, J 

P 

q 

tJ ; 

/ ~ u-iy'- 

-NormaL, 

- Z 

; Z 

. O . -:h 

- Yawing— ,.i.Y 

'Nr 

X >F 

1 

'Yaw i 




S'' 


,s 


Y Absolute coefficients of moment ^ ' 

' n ^ ?> ■ W Af ^ ' -N 

^ ,(roDiiig) (pitchmg) " Cyawmg) 
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41 . PROPELLER SYltf BOLS i 

S P Power, absolute cdejBBcient Cp=- 


p 


' i-, < ' -'I 

u Diameter 

, p Geometric pitcK ■ ^ ' \ - . -rower, aosoiure coemcient Op=— 5 ^ 

/7^ ~ c ■J'** j. t 1 • ' ' y ' j ■ ~ \ ^ ‘ pTv 

piU Pitch ratio - 

y ^Inflow velociity S ' ; , Speed-^power coefficieht=’.W.^Tj 

>, ’■ _ Pi '^'Slipstream^elocityr- S ' ,S \ _/ y -''^Efficiency " ''S', ' ' ' 

f T Thrust, a^olute coefficient > Kevolutions’per second, rps , 

- ' « . TorqS, ibsolWe codScUnt Oo=--fc ' Effeo6T&.helfa angle=tan7'(2^) , , 


-i' Y~ ■' , 'S. numerical RELATIONS. 

1 lb=0.4536 kg . 

- , 1: kg=2.2046 lb 

1 mph=0.4470 mps ' . , 'f , V „ 1 mi= 1,609.35 m= 5^280 ft 

, >1 m==3.2808ft , . ' 


php=76.04 kg-m/s=550 fUb/sed ' 
1 .metric horsepower = 0.9863 bp ^ - 


, JL mps=2.2J69 n^h 
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